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In the standard model (SM), CP violation (CPV) 
arises from the complex phase of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark- mixing matrix Measure- 
ments of the CPV asymmetries in the K and B meson 
systems are consistent with expectations based on the 
SM and, together with theoretical inputs, lead to the de- 
termination of the parameters of the CKM matrix. CPV 
has not yet been observed in the charm sector, where the 
theoretical predictions based on the SM for CPV asym- 
metries are at the level of 10 -3 or below Q. 



In this Letter we present a search for CPV in the decay 
— > K^n^ by measuring the CPV parameter Acp 



defined as: 

_ r(£>+ -> k° s tt+) - r(D- -> jrgTT-) 

CP T(D+ -> K°tt+) + T(D- -> K° s ir-) ' 1 ' 

where T is the partial decay width for this decay. This de- 
cay mode has been chosen because of its clean experimen- 
tal signature. Although direct CP violation due to inter- 
ference between Cabibbo-allowed and doubly Cabibbo- 
suppressed amplitudes is predicted to be negligible within 
the SM [H, K° — K° mixing induces a time-integrated 
CP violating asymmetry of (-0.332 ± 0.006) % 0. Con- 
tributions from non-SM processes may reduce the value 
of the measured Acp or enhance it up to the level of 
one percent [U Q . Therefore, a significant deviation of 
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the Acp measurement from pure K° — K° mixing effects 
would be evidence for the presence of new physics beyond 
the SM. Due to the smallness of the expected value, this 
measurement requires a large data sample and precise 
control of the systematic uncertainties. Previous mea- 
surements of Acp have been reported by the CLEO-c 
((-0.6 ± l.O(stat) ± 0.3(syst))% H) and Belle collabora- 
tions ((-0.71±0.19(stat)±0.20(syst))% 0). 

The data used in this analysis were recorded at or near 
the T(4S) resonance by the BABAR detector at the PEP-II 
storage rings. The BABAR detector is described in detail 
elsewhere [8[. The data sample corresponds to an inte- 
grated luminosity of 469 fb _1 . To avoid any bias from 
adapting the analysis procedure to the data, we perform 
a "blind" analysis where all aspects of the analysis, in- 
cluding the statistical and systematic uncertainties, are 
validated with data and Monte Carlo (MC) simulation 
based on GEANT4 [l(| before looking at the value of 
Aqp- The MC samples include e + e~ — » qq (q = u, d, s, c) 
events, simulated with JETSET [llj and BB decays sim- 
ulated with the EvtGen generator [l2j]. The coordinate 
system defined in [§] is assumed throughout the Letter. 

We select D ± — > K^n^ decays by combining a K® 
candidate reconstructed in the decay mode — > 7r + 7r~ 
with a charged pion candidate. A K® candidate is re- 
constructed from two oppositely charged tracks with an 
invariant mass within ±10 MeV/c 2 of the nominal K® 
mass [ij], which is equivalent to slightly more than ±2.5 a 
in the measured mass resolution. The \ 2 probability 
of the 7r + 7r~ vertex fit must be greater than 0.1%. To 
reduce combinatorial background, we require the mea- 
sured flight length of the K° candidate to be greater 
than 3 times its uncertainty. A reconstructed charged 
track that has pp > 400 MeV/c is selected as a pion 
candidate, where pp is the magnitude of the momen- 
tum in the plane perpendicular to the z axis. At BABAR, 
charged hadron identification is achieved through mea- 
surements of ionization energy loss in the tracking sys- 
tem and the Cherenkov angle obtained from a detector of 
internally reflected Cherenkov light. A CsI(Tl) electro- 
magnetic calorimeter provides photon detection, electron 
identification, and neutral pion reconstruction [8] . In our 
measurement, the pion candidate is required not to be 
identified as a kaon, a proton, or an electron. These 
selection criteria for the pion candidate are very effec- 
tive in reducing the charge asymmetry from track recon- 
struction and identification, as inferred from studying the 
large control sample described later. A kinematic vertex 
fit to the whole decaytree is then performed with no ad- 
ditional constraints [9J. We retain only D ± candidates 
having a % 2 probability for this fit greater than 0.1% and 
an invariant mass m(K®Tr ) within ±65MeV/c 2 of the 
nominal D + mass [3], which is equivalent to more than 
±8 a in the measured D ± mass resolution. Motivated by 
Monte Carlo simulation studies, we further require the 
magnitude of the D ± candidate momentum in the e + e~ 



center-of-mass (CM) system, p*(D ± ), to be between 2 
and 5 GeV/c. This criterion reduces the combinatorial 
background to an acceptable level, but also keeps some 
Z? ± mesons from B mesons decays (they are 8% of 
the selected sample) Additional background rejec- 
tion is obtained by requiring that the impact parameter 
of the Z) ± candidate with respect to the beam-spot 
projected onto the plane perpendicular to the z axis, be 
less than 0.3 cm and the lifetime r 2 ; y (D ± ) be be- 
tween — 12.5 and 31.3 ps. The lifetime is measured using 
L xy {D ± ), defined as the distance of the decay vertex 
from the beam-spot projected onto the plane perpendic- 
ular to the z axis. 

To further improve the search sensitivity, a Boosted 
Decision Tree (BDT) algorithm [14| is constructed from 
seven discriminating variables for each D ± candidate: 
Txy{D ± ), L xy (D^), the CM momentum magnitude 
p*(D ± ), the momentum magnitudes and transverse com- 
ponents with respect to the beam axis for both the 
and pion candidates. Because all the input variables con- 
tains no charge information, no charge bias is expected 
to be introduced by the algorithm and this assumption 
has been verified using a large sample of MC simulated 
events. The final selection criteria are based on the BDT 
output and optimized using truth-matched signal and 
background candidates from the MC sample. For the 
optimization, we maximize the S/^/S + B ratio, where 
S and B are the numbers of signal and background can- 
didates whose invariant mass is within ±31 MeV/c 2 of the 
nominal mass. 

A binned maximum likelihood (ML) fit to the 
TO(i4Tg7r ± ) distribution for the retained Z? ± candidates 
is used to extract the signal yield. The total probability 
density function (PDF) is the sum of signal and back- 
ground components. The signal PDF is modeled as a sum 
of three Gaussian functions, the first two of them with 
common mean. The background PDF is taken as a sum 
of two components: a background from Df — > K®K ± , 
where the K is misidentified as n^, and a combina- 
torial background from other sources. Based on MC 
studies, the yield of D ± — > 7r ± 7r =F 7r ± decays in the fi- 
nal data sample is estimated to be 0.02% of the signal 
and the estimated Acp for this source to be less than 
0.002%. Therefore a PDF to model this component is 
not included in the fit. The background from the decay 
Df K^K* is modeled using a PDF sampled from the 
MC histogram for this mode. The combinatorial back- 
ground is described as a second-order polynomial. The 
fit to the m(i4T!?7r ± ) distribution yields (807 ± 1) x 10 3 
signal events. The data and the fit are shown in Fig. Q] 
All of the fit parameters are extracted from the fit to the 
data sample apart from the normalization of the back- 
ground due to Df — > K^,K ± 1 which is fixed to the value 
predicted by the MC simulation. 

We determine Aqp by measuring the signal yield asym- 
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FIG. 1: Invariant mass distribution for K^-n^ candidates in 
the data (black points). The solid curve shows the fit to the 
data. The dashed line is the sum of all backgrounds, while 
the dotted line is combinatorial background only. The vertical 
scale of the plot is logarithmic. 

metry A defined as: 



where N D +(N D ) is the number of fitted D + — » 
K°ir + (D~ — > K®n~) decays. The quantity A is the re- 
sult of two other contributions in addition to Acp. There 
is a physics component due to the forward-backward 
(FB) asymmetry (App) in e + e~ — > cc, arising from 
7*-Z° interference and high order QED processes in 
e + e~ — > cc. This asymmetry will create a difference in 
the number of reconstructed D + and D~ decays due to 
the FB detection asymmetries arising from the boost of 
the CM system relative to the laboratory frame. There 
is also a detector-induced component due to the differ- 
ence in the reconstruction efficiencies of D + — > K®tt + 
and D~ — > K®ir~ generated by differences in the track 
reconstruction and identification efficiencies for 7r + and 
7r~. While A FB is measured together with A C p using the 
selected dataset, we correct the dataset itself for the re- 
construction and identification effects using control data 
sets. 

In this analysis we have developed a data-driven 
method to determine the charge asymmetry in track re- 
construction as a function of the magnitude of the track 
momentum and its polar angle. Since B mesons are pro- 
duced in the process e + e~ — > T(45) — > BB nearly at rest 
in the CM frame and decay isotropically in the B rest 
frame, these events provide a very large control sample 
essentially free of any physics-induced charge asymme- 
try. However, data recorded at the T(45) resonance also 
include continuum production e + e~ — > qq (q — it, d, s, c), 
where there is a non-negligible FB asymmetry due to 



the interference between the single virtual photon pro- 
cess and other production processes, as described above. 
The continuum contribution is estimated using the off- 
resonance data rescaled to the same luminosity as the 
on-resonance data sample. Subtracting the number of 
reconstructed tracks in the rescaled off-resonance sam- 
ple from the number of tracks in the on-resonance one, 
we obtain the number of tracks corresponding to the B 
meson decays only. Therefore, the relative detection and 
identification efficiencies of the positively and negatively 
charged particles for given selection criteria can be de- 
termined using the numbers of positively and negatively 
reconstructed tracks directly from data. 

Using samples of 8.5 fb _1 on-resonance and 9.5 fb _1 
off-resonance data, applying the same charged pion track 
selection criteria used in the reconstruction of D — > 
K^tt^ decays, and subtracting the off-resonance sam- 
ple from the on-resonance sample, we obtain a sample 
of more than 20 million tracks. We use this sample to 
produce a map for the ratio of detection efficiencies for 
7r + and 7r~ as a function of the track-momentum mag- 
nitude and cos 9, where 9 is the polar angle of the track 
in the laboratory frame. The map and associated statis- 
tical errors are shown in Fig. [2j Since the charm meson 
production is azimuthally uniform, the (j) dependence of 
this ratio is found to be very small and uncorrelated with 
momentum magnitude and polar angle. Therefore, the 
ratio of detection efficiencies is averaged over the (f> coor- 
dinate. The statistical uncertainties can be reduced by 
increasing the control sample size, but this would bring 
a negligible reduction in the final systematic error. In 
the fit procedure described below, the D~ yields, in in- 
tervals of pion-momentum and cos 9, are weighted with 
this relative efficiency map to correct for the detection 
efficiency differences between 7r + and 7r~, leaving only 
FB and CP asymmetries. The average correction factor 
for each interval is —0.09%. 

Neglecting the second-order terms that contain the 
product of Aqp and Apb, the resulting asymmetry can 
be expressed simply as the sum of the two. The parame- 
ter Acp is independent of kinematic variables, while App 
is an odd function of cos 9* D , where 9 * D is the polar angle 
of the candidate momentum in the e + e~ CM frame. 
If we compute A(+\ cos6>J,|) for the D^ 1 candidates in a 
positive cos9* D bin and A(—\cos9* D \) for the candidates 
in its negative counterpart, the contribution to the two 
asymmetries from Acp is the same, while the contribu- 
tion from App has the same magnitude but opposite sign. 
Therefore Acp and App can be written as a function of 
| cos#£>| as follows: 
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FIG. 2: Map of the ratio between detection efficiency for 7r + 
and 7r~ (top) plus the corresponding statistical errors (bot- 
tom). The map is produced using the numbers of tt~ and tt + 
tracks in the selected control sample. 



and 



A CP (\cose* D \) 



A(+\cosd* D \) +A(-\cos6* D \) 



(4) 



Furthermore, the small fraction of the D ± signal yields 
produced from B meson decays have zero FB asymme- 
try. As a result, the measured App from the e + e~ — > cc 
production is slightly diluted, but the Aqp value is unaf- 
fected. 

The selected sample is divided into ten subsamples cor- 
responding to ten cos 9* D bins of equal width and a simul- 
taneous binned ML fit is performed on the invariant mass 
distributions of D + and D~ candidates for each subsam- 
plc to extract the signal yield asymmetries. The PDF 
shape that describes the distribution in each subsample 
is the same as that used in the fit to the full sample, 
but the following parameters are allowed to float sepa- 
rately in each subsample: the yields and the asymme- 
tries for signal and combinatorial events, the mean of the 
second and third Gaussians for the signal PDF, and the 
first order coefficient for the polynomial of the combina- 
torial background. The relative fractions corresponding 
to the second Gaussian are allowed to float only for three 
high-statistics subsamples, while they have been fixed to 
zero for other ones in order to have a converged fit. The 
means of the three Gaussians for the signal PDF, the 
width of the first Gaussian, and the second order coeffi- 
cient for the polynomial of the combinatorial background 
are allowed to float, but they have the same values for all 
the subsamples. Therefore, the final fit involves a total 
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FIG. 3: Acp (top) and Afb (bottom) asymmetries for D — > 
Kgn ± candidates as a function of ] cos 6* D | in the data sample. 
The solid line represents the central value of Acp and the 
hatched region is the ±lcr interval, both obtained from a % 2 
minimization assuming no dependence on | cos#Jj|. 



of 78 free parameters. Using the asymmetry measure- 
ments in five positive and in five negative cosO^ bins, 
we obtain five A FB and five A C p values. As A C p does 
not depend upon cos 8* D , we compute a central value of 
this parameter using a y 2 minimization to a constant: 
Aqp — (—0.39 ± 0.13)%, where the error is statistical 
only. The Aqp and App values are shown in Fig. [3l to- 
gether with the central value and ±1 a confidence interval 
for Acp- 

We perform two tests to validate the analysis proce- 
dure. The first involves generating ensembles of toy MC 
experiments and extracting Aqp for each experiment. We 
determine that the fitted value of the Aqp parameter is 
unbiased, and that the fit returns an accurate estimate 
of the statistical uncertainty. The second test involves 
fitting a large number of MC events from the full BABAR 
detector simulation. We measure Aqp from this MC sam- 
ple to be within ±1 a from the generated value of zero. 

The primary sources of systematic uncertainty are the 
contamination in the composition of particles for the data 
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control sample used to determine the charge asymmetry 
in track reconstruction efficiencies and statistical uncer- 
tainties in the detection efficiency ratios used to weight 
the D~ yields. The charged pion sample selected to de- 
termine the ratio of detection efficiencies for tt~ and n + 
contains a contamination of kaons, electrons, muons, and 
protons at the percent level due to particle misidentifi- 
cation and inefficiencies. This contamination introduces 
a small bias in the Acp measurement due to the slightly 
different particle identification efficiencies between pos- 
itively and negatively charged non-pion particles. The 
particle identification efficiencies, measured in the data 
for positively and negatively charged tracks using the 
method described in the previous paragraphs, are found 
to be in a good agreement with the MC simulation. We 
therefore study this bias using the MC simulated events 
and determine the bias to be +0.05 %. As a result, we 
shift the measured Acp by —0.05 % to correct for the 
bias and then, conservatively, include the same value as 
a contribution to the systematic uncertainty. Therefore 
the bias-corrected value of Acp is (—0.44 ±0.13)%. 

The technique used here to remove the charge asym- 
metry from detector-induced effects produces a small 
systematic uncertainty in the measurement of Acp due 
to the statistical error in the relative efficiency map 
(±0.06%). Using MC simulation, we evaluate an addi- 
tional systematic uncertainty of ±0.01% due to a possi- 
ble charge asymmetry present in the control sample be- 
fore applying the selection criteria. Combining these two 
contributions with the systematic contribution from the 
difference in the composition of the control sample com- 
pared to the signal sample (±0.05%), as described ear- 
lier, the total contribution from the correction technique 
is ±0.08%, which is the dominant source of systematic er- 
ror. We also consider a possible systematic uncertainty 
due to the regeneration of K° and K° mesons in the ma- 
terial of the detector. K° and K° mesons produced in the 
decay process can interact with the material around the 
interaction point before they decay. Following a method 
similar to that described in [l^j], we compute the proba- 
bility for K° and K° to interact inside the BABAR track- 
ing system. We numerically integrate the interaction 
probability distribution, which depends on the measured 
nuclear cross-section for (assuming isospin symme- 
try) , the amount of material in the BABAR beam-pipe and 
tracking detectors, the K°/K° time evolutions, and the 
K® kinematic distribution and reconstruction efficiency 
as determined from simulation studies. From the dif- 
ference between the interaction probabilities for K° and 
K°, we estimate a systematic uncertainty of ±0.06%. Mi- 
nor systematic uncertainties from the simultaneous ML 
fit are also considered: the choice of the signal and back- 
ground PDF, the limited MC data set to estimate the 
normalization of Df — > K^K^, and the choice of bin- 
ning in cos#|), for a total contribution of ±0.01%. The 
combined systematic uncertainty in the CP asymmetry 



measurement including all the contributions is calculated 
as the quadrature sum and is found to be ±0.10%. 

In conclusion, we measure the direct CP asymmetry, 
Acp, in the — > K^n^ decay using approximately 
800,000 D ± signal candidates. We obtain 



A C p = (-0.44 ± 0.13 ± 0.10) %, 



(5) 



where the first error is statistical and the second is sys- 
tematic. The result is consistent with the prediction of 
(-0.332 ± 0.006)% for this mode based on the SM. 
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